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Abstract 



D , We compute critical exponents in a Z2 symmetric scalar field theory in three dimensions, using 



Wilson's exact renormalization group equations expanded in powers of derivatives. A nontrivial 
relation between these exponents is confirmed explicitly at the first two orders in the derivative 



H ' expansion. At leading order all our results are cutoff independent, while at next-to-leading order 

they are not, and the determination of critical exponents becomes ambiguous. We discuss the 
possible ways in which this scheme ambiguity might be resolved. 
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The exact renormalization group [^Q is in many instances a powerful analytical tech- 
nique for the study of nonperturbative renormalization group flows. In practice one must 
always truncate an infinite tower of equations, and in the past surprisingly good results for 
critical exponents and renormalization group flows of scalar fleld theories m2 < d < 4 have 
been obtained [l5|-|To| by suppressing all the equations except that for the effective potential. 



There have been several attempts |jTT|-|T^ to improve the accuracy of these calculations by 
also keeping terms in the effective action containing two derivatives of the flelds, and thus 
allowing for wave function renormalization. In this letter we will consider the important 
issue of scheme dependence within this context. 

A Euclidean field theory may in general be regularized by employing a certain cutoff 
function Ka{p), which suppresses modes with momenta higher than some scale A. The 
exact renormalization group equations are then constructed so as to ensure that as A is 
reduced, and thus more and more modes are integrated out, the vertex functions of the 
theory are varied in such a way that the Green's functions remain unchanged. This provides 
a nonperturbative definition of the theory. Now in principle physical observables should 
be determined unambiguously for a broad class of cutoff functions Ka{p)- This is what we 
mean by scheme independence. Indeed in certain truncation schemes this independence can 
be proven; for instance, it can be shown that S-matrix elements are scheme independent 



to any given order in perturbation theory [|T3| , even though beta-functions and anomalous 



dimensions are in general scheme dependent beyond leading order. However in other 
truncations scheme independence cannot be taken for granted. In the derivative expansion, 
in particular, this problem has not so far been considered: a particular choice of cutoff 
function is made (e.g., sharp cutoffs [^0, power-like cutoffs [jl^, or exponential cutoffs 



T^,|I^), giving perhaps the impression that the results obtained are in some sense unique. 

Here we will compute various critical exponents for a Z2-symmetric scalar field theory 

in (i = 3 dimensions using the simple and elegant exact renormalization group equations due 

to Wilson [jl| and Polchinski [|| : previous authors [|5|-p|, p!3| , p^ have based their computations 



on the Wegner-Houghton and Weinberg equations. We work to all orders in the fields, 
but expand the (Wilson) effective action in powers of derivatives. We will show that all the 
scheme dependence can be absorbed into 2n parameters at n-th order in the expansion. 
At leading order, n = 0, the critical exponents are then scheme independent: they are 
independent of the particular choice of cutoff function. However at next-to-leading order 
the results become scheme dependent. This makes it particularly difficult to compute 
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unambiguously the critical exponent ry (corresponding to wave function renormalization) . 
We will attempt to control this dependence by employing a minimal sensitivity criterion. 
We also discuss the 'scaling' relation between the magnetic deformation exponent Xh 
and 77, showing explicitly that this is satisfied exactly (in any scheme) at the first two 
orders in the derivative expansion. 

1. The Exact Renormalization Group 

Exact renormalization group equations are normally obtained in two steps: incomplete 
integration of modes, followed by rescaling. The first step is most easily described in the 
path integral approach developed by Polchinski Q. Consider the following Wilson action 
for a scalar field theory 

S[^;A] = i [^p^-pP^\p'')+S,^t[^;Al (1.1) 

Jp 

where J = J (2-J)d ' -Pa(p^) is an analytic function [T^] with a single pole on the negative 
real axis of unit residue, which falls sufficiently fast as p^/A^ -^ 00 that all modes with 
Euclidean momenta p^ ^ A^ are suppressed, and Si^t contains all higher order interactions. 
We will write generically, 

P,(,^) . im^ , (1,) 

where K{p^ / K'^) is the cutoff function responsible for damping high momentum modes. 
The exact renormalization group then makes use of two momentum scales, an ultraviolet 
regulator scale Aq, which is held fixed, and A which controls the scale down to which high 
momentum modes have been (loosely speaking) integrated out. a What is of relevance for 
the flow is their ratio or, equivalently, t = In -^; renormalization group flows then go from 
the far ultraviolet (t = 0, A = Aq) to the far infrared (t -^ 00, A ^ 0). As A is reduced from 
Aq to zero the interaction, 5'int [4>] A] , must be evolved in such a way that the amputated 
Green's functions of the theory remain unaltered; this can be achieved if 5'int satisfles the 
renormalization group equation 

.rf„ 1 /■.<WW) /■<ssi„,<ssi„, .s"Si„t ^ 1 „ /■ p_i, 2, ,, „ 



More properly we should also include a mass scale m <C Aq (the position of the zero in Pa); 
while important in the extreme infrared (A <C m), it may be ignored in the vicinity of the Wilson 
fixed point, since there m/A will be very small. 



(Contact terms have been suppressed here; see ref.[15] for more details.) Here r]{t) is the 
anomalous dimension of the (renormalized) field (fpi 

A^<^p = -i^Wv'p- (1-4) 

The above renormalization group equation (|1.3|) may be further simplified by rewriting 
it in terms of the complete Wilson action S[f; A]: suppressing the contact term this gives 

dA "^ Jp dA \5(pp5(p-p 5(ppScp-p ^ ^ScppJ' 

The explicit 77 dependence has now dropped out because renormalizations of the fields 
can be absorbed directly into corresponding renormalizations of the vertices. As we shall 
see, this equation []I| is rather easier to use than the other exact renormalization group 
equations [^^ employed in p|-p|, p!3| , p^ , due essentially to the relative simplicity of its non- 
linear term. 

Finally we consider the rescaling step. For this purpose we have to write the cou- 
plings in 5" in a manifestly dimensionless way, since it is dimensionless couplings which 
parametrize the theory and whose flow we should study. Thus, we write 

-^int = Yl / s^{pi,...,pr,;t) ipp^---0p^ 5{pi + ...+Pn) , (1.6) 

n=even "^Pi-'-P" 

where pi = A~^pi, cpp = A^+^'^c/jp, and the couplings (or 'vertex functions') 
SniPiT ■ ■ ■ ,Pn]i) are dimensionless. We can then rewrite ( p..3|) as a set of equations for 
the flow of the couplings Sn- 



S=l K'ip^){jl^--J^]+dS 



d m o^2K'{p')\,, SS f .^ ^ d' 5S 



^i\'-i-^-'^'im)'^Wp-J/^' ops,.. 



;i.7) 



where S is given by a similar expression to ( |1.6| ), but with the dimensionless couplings Sn 
replaced by their partial derivatives with respect to t, and the prime in the momentum 
derivative means that it does not act on the momentum conserving delta functions of ( |1.6| ) . 
The form of r]{t) along the flow is not yet determined: it will depend on the chosen 
normalization of the fleld (fp. For the bare (unrenormalized) flelds used in 0,|15[ r]{t) = 0, 
and the interaction S'lnt will develop quadratic terms which change the residue of the pole in 



the propagator. Alternatively one can choose the anomalous dimension -q^t) in such a way 
that no such terms arise, and the residue thus remains fixed, by imposing the normalization 
condition 



, d d 5'^S 
A—T- 



dA \ dp^ 5(/?p(/?_ 



0, (1.^ 






at each point of the renormalization group flow. This then implicitly deflnes the anomalous 
dimension rjit), which at flxed points of the flow becomes equal to the critical exponent 

2. Derivative Expansion and Scheme Dependence 

The overwhelming complexity of the full exact renormalization group equations means 
that in practice some sort of truncation must be employed. Since the action may be 
naturally expanded in powers of the flelds and their derivatives, it is natural to consider 
truncations of either or both of these expansions. Truncations in the number of powers 
of the flelds, as employed in j^^] seem to be poorly convergent, with many spurious 
solutions (see |T^ for a detailed analysis), although the convergence seems to improve if 



the expansion is made about the minimum of the potential P,|l^ • A better procedure [l5|,[Tl|- 



|T^ seems to be to work to all orders in the number of flelds, and expand the interaction 
only in powers of derivatives. In this way it is hoped that the essential features of the 
long-wavelength physics will be retained. 

Truncating at second order in derivatives, we thus write 

^int = A'^ /" A {v{ip,t) + z{ip,t)K-\d^ipf + • ■ •) (2.1) 

or, alternatively, 

with v{ip,t) = '^Vn{t)0^, z{ip^t) = J2^n{i)0^i and where now zo{t) = due to the 
normalization condition ( p..8|) . The exact renormalization group equation (|1 . 7|) can now 



be projected onto a series of equations, order by order in p , in which (f is taken to 



be effectively a constant (and denoted by x). This procedure is somewhat lengthy but 
unambiguous, and results in an infinite set of coupled partial differential equations:0 

i, = Iqv" + 2hz - Kqv'^ + (1 - f - ^)xv' + dv , 

z = loz" + Kiv"^ - 4Kozv" - 2Koz'v' + (1 - f - f )x/ - r?^ - f , (2.3) 



where dots and primes denote partial derivatives with respect to t and x respectively. The 
number of equations grows rather fast with increasing powers of p'^ ; there are three more 
equations at the next order. So far, there seems to have been no attempt to go to order 
p'^; here we will just study the two equations ( |2.3D . The numbers K^, 1^, n = 0, 1, . . . are 
remnants of the cutoff function K{z) and thus parametrize the scheme dependence of the 
equations: 

Kr, = (-)"+1k("+i)(0), Ir, = - [{pTK'ip^) = -Qd rdzz''/^-^+^K'{z) , (2.4) 

Jp Jo 

where K^'^'^ is the n-th derivative of K, and fid = 2/(r(|)(47r)'^). At zeroth order in 
p^, all the cutoff dependence is parametrized in terms of Kq and /q, while at first order 
Ki and /i are also needed; at n-th order we need Kn and In- Thus the further we go 
in the derivative expansion, the more information about the cutoff function is needed; 
short-distance properties of it only gain relevance as the derivative expansion is carried to 
higher orders. Incidentally, we note that in a calculation of perturbative renormalization 
group functions in this setting, these same remnants of the cutoff function are present and 



determine the scheme dependence of some of the results [^ | . 

We have just shown that the cutoff dependence of eq. (|2.3|) is given in terms of at 
most 2?T, + 2 arbitrary parameters at order p'^^. In fact we can exploit global rescalings of 
the variable x and of the function v to further simplify the first two equations to 

v = v" + 2Az -v'^ + (l-^- ^)xv' + dv , 

(2.5) 
z = z" + Bv"^ - Azv" - 2z'v' + (1 - f - l)xz' - vz - ^ , 



where 



, hKo Ki , , 



The zeroth order equation (with z = ?] = 0) seems to have been first derived in [p^ ; different 



versions of the first order equations may be found in [11,12|. 



It follows that at zeroth order {z = 0) there is no cutoff dependence at all, while all 
cutoff dependence at order p'^ is reduced to a two-parameter family {A, B). Furthermore, 
reasonable cutoff functions (for example, exponentials of polynomials [|T5| ) will decrease 
smoothly and monotonically, so all moments In will be positive and of similar magnitude; 
the derivatives at the origin, Kn, will also be loosely of the same order. In other words, A 
and B will both be positive and loosely of order unity. For example, with an exponential 
cutoff Oig] K{z) = e""^, A = d/2 and B = 1. On the other hand a sharp cutoff 
K{z) = 9(1 — z), corresponds to a singular point in the A-B parameter space, since A = 
while B is undetermined, and can thus only be made well-defined through some limiting 
procedure (as discussed in [§). The power-like cut-oG K{z) = (1 + z'^)~^ employed in [^| 
is also singular since A = and B is infinite. In these cases it becomes necessary to use 
other more complicated exact renormalization group equations [|l-0. 

3. Leading Order Numerical Results 

The truncation of the derivative expansion of the exact renormalization group equation 
at zeroth order is obtained trivially from eq. ( |2.5| ) by writing z{x, t) = and r] = 0. Setting 
r] to zero is justifed only by the a priori suppression of all momentum dependent terms in 
the interaction; it turns out to be a good approximation at zeroth order because the actual 
value of r]^ is very small. Writing f = v' (to eliminate the arbitrary field independent part 
of v) the leading order exact renormalization group equation is then 

/ = r - 2//' + (1 + f )/ + (1 - l)xf' . (3.1) 

This is considerably simpler than equivalent truncations of other exact renormalization 
group equations 0,|^, since all non-linearity is reduced to a quadratic term which is easily 
controlled numerically. Furthermore, eq. ( p.l|) is manifestly scheme independent (that is, 
independent of the particular choice of cutoff function). This is of great importance since 
the results we will display shortly are the leading orders of the expansion. 
We begin by solving the fixed point equation 

= //' - 2fj: + (1 + f )/, + (1 - l)xfj , (3.2) 

with initial conditions /*(0) = (by the Z2 symmetry) and /*'(0) = 7. We now concentrate 
for definiteness on dimension d = 3 where a comparison with other results is readily 

6 



available. The nontrivial solution at the ultraviolet fixed point only exists for one value of 
the initial condition 7 = 7*. Numerically, this is found by tuning 7 so as to let /* run as 
far as possible; in this way we determine 7=^ = —.228601293102 . . . (where all figures are 
significant: the behaviour of f*{x) at large x is extremely sensitive to the precise value of 
7*, always diverging at some finite value of x for 7 7^ 7^,). 

It is in fact trivial to determine the asymptotic behavior of the solution f^{x) to 
( P?2|) at large x: /*(x) ~ x, and thus v^{x) ~ x^. This is in marked contrast to the x^ 
asymptotic behavior of solutions to the equations considered in P JT^] , resulting from their 
rescaling terms. This difference stems from the different structure of the non-linearities in 
the equations; we find it intriguing since it is ultimately responsible for the differences in 
the numerical results. 

Departures from the fixed point solution are governed by operators whose critical 
exponents are universal. In our case, the eigenvalue equation to be solved is given by the 
linearization of eq. (|3J]): writing /(x, t) = f*{x) + J2n '^nfi'n(a;)e^"*, with a^ small, 

g'; - 2if,g'^ + fjg^) + ^g^ - y^ = X^g^ , (3.3) 

where gn is the n-th eigenfunction with eigenvalue Xn- The critical exponent of the only 
relevant, Z2-symmetric, operator at the infrared fixed point is known as z/ = 1/Ai, while the 
first irrelevant one is known as cu = — A2. We investigate eq. (|3.3|) numerically and find the 
results given in Table III below (LO). Further eigenvalues for irrelevant operators may be 
determined similarly, but become increasingly meaningless as they depend more and more 
on short-distance dynamics. Curiously our leading order result u = 0.649 is rather closer 



to the values obtained using more traditional methods [^ than that obtained using the 
leading order approximation to other exact renormalization group equations (i.e. u = 0.690 
for a Wegner- Houghton equation [|^-0], u = 0.660 for an equation for one particle irreducible 
vertices with a power cut-off |]1^). However before we can attach any significance to this 
observation we need to determine the size of the error induced by the truncation of the 
derivative expansion, and this we can only do by calculating the next order correction. 



4. Next-to-Leading Order 

The main drawback of the simple zeroth order approximation described above is the 
absence of wave function renormalization: 77 is set to zero. Equation ( p.l|) thus includes 
only the physics of zero modes, whereas wave function renormalization is precisely related 



A?7 X 10^ 


B = 0.25 


0.5 


0.75 


A = 0.5 


-0.76 


- 2.8 


-6.1 


0.7 


-0.28 


- 0.9 


-1.9 


0.8 


-0.03 


0.0 


0.2 


0.9 


0.21 


1.0 


2.5 


1.1 


0.71 


3.0 


7.1 



Table I: The difference Arj = r]2 — rji for a range of values of the scheme 
parameters A and B. 



to the kinetic term which carries derivatives. This leads us to consider the next-to-leading 
corrections, of order p'^ in the derivative expansion. 

To study numerically the order p^ approximation of the exact renormalization group 
equation eqns. ( J2.5D we first consider the fixed point differential equations /(t, x) = 

0, i(t,x) = 0: 

= /:' + 2A< - 2fj: + A+/, + A-xfi , 

(4.1) 
= z':+ Bfi^ - Azfi - 2z'f, + A-xz', - r^z, - f , 

where A^ = 1 ± d/2 — ri/2. This pair of non-linear coupled ordinary differential equations, 
together with appropriate initial conditions, determine the shape of the Wilson fixed point 
action at order p^. Three of the initial conditions are already fixed: two by invoking Z2 
symmetry, /h<(0) = and z'^{Q) = 0, and the third by the normalization condition eq. ( |1.8|) , 
z^{Q) = 0. For any fixed values of the parameters A and S, a well-behaved solution of eqs. 
( p75| ) {i.e., a solution defined for all finite x) only exists for the unique value ry = 77* of the 
anomalous dimension and a unique value 7* of the initial condition /^(O) = 7. For any 
other choice of 77 and /^(O) the solutions /*(x) and z^{x) end at singularities at finite x. 
By contrast, the correct solutions behave asymptotically, as x — > cxd, as 



/,(x)~(l-77,/2) x + A;x-^^/^ 
z^[x) ^ — ^ ^ -f k 



2k' 



A+(2-r7, 
2) - 77, 



X 



-2A+/A- 



■^ + 0(x 



-3A+/A~-2^ 



v* 



A+ 



X 



-A+/A- 



v* 



-^ + 0{x 



-2A+/A- 



(4.2) 
where fc is a constant. 

Equations ( fl.l|) are nonlinear and stiff. Furthermore, finding the correct solution 
involves a double fine tuning in 77 and /^(O). Thus it becomes extremely difficult to simply 



integrate up the equations as we did in the zeroth order case, and it is much easier to 
proceed recursively. So first we set z^{x) and rj to zero and find the unique value 70 of 
/'(O) for which the first equation ( p] ) has a non-singular solution, /o(x), just as we did 
for the zeroth order approximation. Now we use this solution in the second equation ( |4.1J ) 
and solve for zi{x) by fine-tuning rj to the appropriate value rji. With the values of 771 and 
zi{x) we return to the first equation and obtain a new solution /i(x), with /{(O) = 71, 
and so on. In this way we obtain a sequence of functions /o, -Zi, /i, ,22, /2, • • ■ and constants 
705 ^i5 7i5 ^2, 72, • • • which (at least for a reasonable range of A and B) converge rather 
rapidly to the exact solution /h<(x), ^^(a;), 7*, ^*- Solutions /*(x) and z^{x) are shown in 

fig- I- 

Again for d = 3, a numerical analysis of this sequence of solutions shows that the 
rate of convergence of the iterative procedure is in effect controlled by A. Indeed for any 
fixed B, the convergence is fastest (in the sense that corrections to r] in sucessive iterations 
are minimized) when A = 0.8. Moreover, the function r]{A) has an infiexion point at 
this particular value of A. This is demonstrated in Table I where we show the difference 
Ar] = r]2 — r]i for different values of the parameters A and B. We will use this criterion to 
fix, from now on, the value of A to 0.8. 

Within the accuracy of our computations, we could not find a similar minimal sen- 
sitivity criterion for the parameter B either in the fixed point solutions or in the critical 
exponents. The values of r]^ are given in Table II. The dependence of the anomalous dimen- 
sion 77* in the parameter B is almost linear, as can be inferred by the following argument. 
By varying the initial conditions of /* and 2;* at x = in the first equation ( ^.5| ) and 
assuming a (fast) convergence of the iterative process, we can relate the second derivative 
of Zi at the origin to the parameter 79. Using this relation and the projection onto a; = 
of the second equation ( |2.5D one can deduce the approximate relation: 

7*' 
V* ^ 2S-^ , (4.3) 

in good qualitative agreement with our numerical results. Although no serious results can 
be obtained from this estimate of 77*, it is interesting that it has no stationary points in 
either A or B. The optimal value oi A = 0.8 is, thus, due to higher order iterations. 
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V 


V 


LO 


B = 0.25 


0.01898 


0.63706 


0.7016 


0.5 


0.03762 


0.62522 


0.7733 


0.75 


0.05595 


0.61595 


0.8501 



Table II: Critical exponents at next-to-leading order for A = 0.8 and a range 
of values of B. 

However, to the precision we work with, the dependence of r]^ on B remains almost linear 
even when subsequent iterations are considered. 

Once a fixed point solution is found we can study the critical exponents associated 
with this solution, by again linearizing equations ( p.5|) around the fixed point. This gives 
rise to the eigenvalue problem 

9nix) - 2g^{x)fJ{x) - 2f^{x)g'^{x) + A'^g^{x) + A~xg'^{x) + 2Ah'^{x) = Xn9nix), 

Ki^) - '^f*ix)h'^{x) - 2gnix)z/{x) - 4f/{x)hn{x) - 4:g'^{x)z^{x) 

+ A~xh'^{x) - Tj^hnix) + 2BfJ{x)g'^{x) = X^h^ix), 

(4.4) 

where /*(x) and z^{x) are the fixed point solutions, 77* the corresponding anomalous dimen- 
sion, and {gn{x) , hn{x)) is the eigenvector (scaling operator) corresponding to the critical 
exponent A^. 

For Z'2-symmetric perturbations, we now use as initial conditions gn{0) = 0, h'^{0) = 0, 
which are fixed by the Z2-symmetry, and fi'^(O) = 1 as a normalization. Solutions are again 
found by an iterative procedure, where the two parameters to be tuned are now An and 
hn{0) = hoB When (i = 3, we find as expected that for A = 0.8 and a reasonable range 
of values of B all the eigenvalues were very close to those found at zeroth order. Results 
after four iterations, for various values of B, are given in Table II. 



5. The Magnetic Deformation Exponent 

Thus far we have only considered Z2-symmetric perturbations about the fixed point. 
It is also interesting to consider Z2-antisymmetric deviations from the fixed point, since 



Although ho is very small, it is nonzero due to the fact that we have implicitly relaxed the 
normalization condition eq. ( |1.8D by setting r]{t) = r]* away from the fixed point. 
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V 


V 


ijj 


LO 





0.649 


0.66 


NLO 


0.019 - 0.056 


0.616 - 0.637 


0.70 - 0.85 


NLOS 


0.035 


0.627 


0.76 


Ref.|l| 


0.030 - 0.040 


0.630 - 0.631 


0.75 - 0.85 



Table III: Critical exponents at leading order (LO), next-to- leading order 
for B in the range .25 - .75 (NLO) , next-to-leading order with the scheme 
dependence fixed to give 77 correctly (NLOS), and finally the average value 
obtained from other methods [19|. 



these yield another critical exponent Xn- In the analogous Ising system, such perturbations 
measure the response of the system to a weak external magnetic field. 

To find these antisymmetric eigenvectors to eq. ( [1.4| ), we use as initial conditions 
fi'n(O) = 1,(7^(0) = 0, /in(0) = and /i^(0) = /iq very small. This does not give an 
independent critical exponent for the theory, but it does allow us to verify for (i = 3 the 
'scaling relation' |]n_ 

(5.1) 






giving us a nontrivial consistency check on our results for r\^,. At order jp and at order 
p^, with A = 0.8 and several values of S, we have calculated A^, and have confirmed 



numerically that the relation (|5.1|) is always satisfied. 

Indeed it is possible to prove (|5.1| ) rigorously for any d, any scheme, and at leading or 
next to leading order in the derivative expansion: it is straightforward to verify that 



(^.(x),/i.(x)) = (/:(a:)-(l-^),<(x)) 



(5.2) 



is an exact eigenvector of eq. (| 
Incidentally, 



with eigenvalue A^ given by eq. ( |5.1| ). 



(^„(x),/i„(a;)) = {f^{x),z^{x)) 



(5.3) 



is also an exact eigenvector, with eigenvalue A/ = — 1-|-|-(-|. At(i = 3, there are 
thus three relevant eigenvalues: At = 1/^^, ^h and A/. Just as \u and At correspond 
to perturbations proportional to the operators (/? and (/?^, respectively. A/ corresponds to 
perturbations containing the operator (/?'^, but proportional to the equations of motion. 



11 



6. Conclusions 

Like perturbation theory, the derivative expansion to any finite order represents a 
truncation of the renormalization group. We find it interesting that our results for anoma- 
lous dimensions are scheme independent (or, more precisely, cutoff independent) to first 
nontrivial order and scheme dependent beyond that. 

Our results for the leading critical exponents of a scalar field theory with Z2 symmetry 
in three dimensions at the Wilson fixed point are summarized in Table III, together with 
values obtained from different methods ||19|| (resummed perturbation series, e-expansions. 



lattice simulations, etc.). At leading order our (scheme independent) equation seems to 
do remarkably well; essentially this is because the true value of r\ is very small. At next- 
to-leading order the results improve in the right direction (and are consistent with those 
in 



r3| , [r^ ) , but scheme dependence makes it impossible to give a precise estimate of the 
errors at this order. Our attempts to fix the scheme by some sort of 'minimal sensitivity' 
criterion were only partially successful, because 77 depends approximately linearly on the 
scheme parameter B. A next-to-next-to-leading order computation might give results with 
a weaker scheme dependence, allowing perhaps for a determination of r\ and thus a sensible 
estimation of the errors incurred in other quantities also. However, at present scheme 
dependence severely limits the accuracy of the technique, in particular for quantities which 
depend sensitively on 77. The scheme ambiguity of our results may however be significantly 
reduced if we simply choose a scheme which gives the 'known' value for r\: our results 
for V and u) are then both competitive and consistent with the 'known' results for these 
exponents. 

We also discussed two other relevant eigenvalues, corresponding to asymmetric pertur- 
bations about the fixed point. We showed that the scaling relations relating these relevant 
eigenvalues to the other critical exponents are satisfied exactly by our truncated set of 
equations, and thus give a useful consistency check on the numerical calculations. 

It should be possible in the future to use these truncated exact renormalization group 
equations to compute other quantities of physical interest (and in particular Green's func- 
tions), in a wider range of models. In particular it would be possible [^ to include (chiral) 
fermions, and thus study (for example) the electroweak sector of the standard model (with 
gauge couplings suppressed). 
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Figure Captions 

Fig. 1. The fixed point solution at next to leading order for A = 0.8 and B = 
0.25,0.5,0.75: a) /*(x) b) z^{x). Ttie (sctieme independent) leading order so- 
lution is also shown (dashed). 
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